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ABSTRACT 

Spitzer Space Telescope spectra of the low mass young stellar object (YSO) IRS 46 (Lboi ~ 0.6 L0) in 
Ophiuchus reveal strong vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2. This is the only 
source out of a sample of '^100 YSO's that shows these features and the first time they are seen in the spectrum 
of a solar-mass YSO. Analysis of the Spitzer data combined with Keck L- and M-band spectra gives excitation 
temperatures of > 350 K and abundances of 10"^- 10"^ with respect to H2, orders of magnitude higher than 
those found in cold clouds. In spite of this high abundance, the HCN 7 = 4-3 line is barely detected with the 
James Clerk Maxwell Telescope, indicating a source diameter less than 13 AU. The (sub)millimeter continuum 
emission and the absence of scattered light in near-infrared images limits the mass and temperature of any 
remnant collapse envelope to less than 0.01 Mq and 100 K, respectively. This excludes a hot-core type region 
as found in high-mass YSO's. The most plausible origin of this hot gas rich in organic molecules is in the inner 
(< 6 AU radius) region of the disk around IRS 46, either the disk itself or a disk wind. A nearly edge-on 2-D 
disk model fits the spectral energy distribution (SED) and gives a column of dense warm gas along the line of 
sight that is consistent with the absorption data. These data illustrate the unique potential of high-resolution 
infrared spectroscopy to probe organic chemistry, gas temperatures and kinematics in the planet-forming zones 
close to a young star. 

Subject headings: infrared: ISM - ISM: individual (IRS 46) - ISM: jets and outflows - ISM: molecules - 
planetary system: protoplanetary disks - stars: formation 



1. INTRODUCTION 

The presence of gas-rich disks around young stars is well 
establis hed observationally and theoretically (see review by 
lGreavesii20 05). but comparatively little is known about their 
chemical structure. A good understanding of the chem- 
istry is important since part of the gases and solids in pro- 
toplanetary disks will end up in future solar systems where 
the y may form the basis for prebio ti c species (see reviews 
bv lEhrenfreund & CharnlevI llOOd iMarkwick & CharnlevI 
120041) . Virtually all observational studies of molecules 
other than CO have been limited to the outer regions, 
where simple organic molecules such as HCO^, CN, HCN, 
and H7CO have been detected at millimeter wavelengths 
(e.g. iDutrev . Guilloteau & Guelin"1997*. 'Kastner et al"l99'T; 
lOietalJ 12003: Thi, van Zadelhoff & van Dishoeck 2004). 
Because of beam dilution, these observations cannot probe 
radii < 50 AU from the star, which is the relevant zone for 
planet formation. High-resolution infrared (IR) spectroscopy 
has found CO emission from the warm dense gas in the in- 
ner d i sk region (Naiita. Carr & Mathieu 2003; Brittain et al. 
|2003; Blake & Boogert 2004), but H2O is the only molecule 
besides CO and H2 which has been convincingly detected 
(|Carr, Tokunaga & Naiita 2004). 

Models of inner disk chemistry initially focused on our own 
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primitive solar nebula (see review bv lPrinDlI99 3') but now also 
consider exosolar systems. They have grown considerably 
in sophistication, including non-equilibrium chemistry, gas- 
solid interactions, radial and vertical mixing, and the effects 
of UV radiation and X-rays from the centra l star on the gas 
temperature and molecular abundances (e.g., Markwick et al 
2002; Gail 2002; Ilgneret al. 2004; Glassgold, Naiita & Ige" 
^04; Gorti & Hollenbach 2004). Large abundances of or- 
ganic molecules like C2H2 and HCN are predicted in some of 
these models in the inner few AU, but no observational tests 
have been possible to date. 

The sensitive Infrared Spectrograph (IRS) on board the 
Spitzer Space Telescope opens a new window to study 
molecules in disks through IR pencil-beam line-of-sight ab- 
sorption spectroscopy. The Spitzer c2d legacy pro gram "From 
Molecular Cores to Planet-Forming Disks" ( Evans et"al] 
2003) is collecting a coherent sample of IRS spectra of low- 
mass YSO's in five nearby star-forming regions. To date 
more than 100 Class I and Class II sources have been ob- 
served. Of these only one source, IRS 46, shows strong 
gas-phase absorption bands of hot molecules. These gas- 
phase bands have previously been seen only toward deeply 
embedded high-mass YSO's, where they have been associ- 
ated with the inner (< 1000 AU) warm dense reg ions of the 
spher i cal envelopes, also kno wn as 'hot cores' JCarr et alJ 
[19951 lT,ahuis & van Dishoeckll2?M iRoonman et al.ll2003^ ■ 
The sources studied with Spitzer have factors of 10"*- 10^ 
lower luminosity, thus limiting the maximum temperatures 
and amount of warm gas in the envelope. 

IRS 46, also known as YLW16b and GY274, is part 
of the Ophiuchus cloud at a distance of -^125 pc 
ide Geus. de Zeeuw & LubllI989h . It is classified as a Class 
I source, i.e., a protostar with a compact accretion disk em- 
bedded in a more extended a nd collapsing envelop e, based on 
its near- and mid-IR colors dAndre & Montmeri3ll994l) with 
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Fig. 1 . — Composite SED of IRS 46 from 1.25 fim to 1.3 mm including 
the full Spitzer IRS spectrum from 10-37 ^tm, 2MASS J.H.Ks imaging and 
photometry (Fj = 0.33 mJy, Fh = 9.03 mJy, Fk., = 77.5 mjy), a VLT-ISAAC 
L-band spectrum at R 1200 (Pontonnidan et al. 2003), an ISOCA M-CVF 
5 — 16 fim spectrum atR^ 50, FgsQ = 150 mJy (see Ridge et al. 2006), and 
^1.3 mm = 28 mJy I Andre & Montmerle 1994). Overplotted in grey is a SED 
disk model of a nearly edge-on self-shadowed flaring disk of 60 AU radius 
with L(,oi RJ O.6L0 and a disk mass of ~ O.O3M0 seen at an inclination 
of ~ 75°. The insert shows a zoom-in on the mid-lR region showing the 
observed C2H2, HCN, and CO2 molecular absorption features and the CO2 
ice band at 15 ^m. 
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Fig. 2. — Blow-up of the IRS 46 normalized Spitzer-IRS spectrum cover- 
ing the C2H2 Us = 1-0, HCN U2 = l-Q, and CO2 1^2 = 1 - bending mode 
ro- vibrational absorption bands. Included in grey is a best fit synthetic spec- 
trum. 

Lbol ~ 0.6 Lq jBontemps et alJ2001l) . However, the complete 
SED is also consistent with a Class II source viewed nearly 
edge-on, i.e., a pre-main sequence star with a disk but with- 
out a significant collapse envelope (see Sect.|4}. It is similar 
to the SED of t he nearly edge-on disk CRBR 2422.8-3423 
iPontoppi dan et al. 20Q5.) . but IRS 46 has less dense fore- 
ground material. The favorable inclination of CRBR 2422.8- 
3423 allows for the study of ices in the outer region of the 
circumstellar disk. IRS 46, also profiting from a favorable in- 
clination, may prove to offer a unique opportunity to directly 
study the gas temperature and hot gas-phase organic chem- 
istry in the inner disk. 

2. OBSERVATIONS 

IRS 46 was observed with Spitzer-IRS in the SH (9.9- 
19.6 /im) and LH (18.7-37.2 fim) medium-resolution modes. 
The observation was positioned at 16''27'"29^4 -24°39'16".3 
(J2000) and acquired on 2004 August 29 as part of AOR# 
0009829888. Low resolution SL and LL observations are 
scheduled as part of the Spitzer GTO program but have not 
yet been taken. Data reduction started from the BCD images 
using S 11. 0.2 Spitzer archive data. The processing includes 
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Fig. 3.— Top: Keck-NIRSPEC L-band spectrum showing HCN z/3 C- 
H stretching mode absoiption lines with tickmarks at the rest wavelengths. 
Bottom: HCN and CO velocity components derived from Keck-NIRSPEC L- 
and M-band spectra using absorption lines with an atmospheric transmission 
> 40%. Plus symbols show the observed data, solid lines the fitted Gaussian 
velocity profiles. In the '^CO profile plot the HCN profile is included (dashed 
grey line) showing that the HCN and CO lines are blue-shifted by a similar 
amount but that the CO profile is broader. 

bad-pixel correction, extraction, defringing, and order match- 
ing using t he c2d anal ysis pipeline (Lahuis, in preparation; 
Kessler-Sil acci et alJ20'05t) . 

Figure shows the SED composed from the full Spitzer- 
IRS spectrum and complementary archival and literature data. 
Figure|2]shows the part of the (normalized) IRS spectrum re- 
vealing the C2H2 1^5 = 1-0, HCN v2=l-0, and CO2 V2=l-0 
bending mode ro-vibrational absorption bands. Included is a 
best-fit synthetic spectrum (see Section|3}. 

To constrain the source size of the warm gas, JCMT^ HCN 
7 = 4- 3 and CO 7 = 3- 2 spectra at 354 GHz and 345 GHz 
were obtained using receiver B3. A weak HCN line with 
Tmb = 0.035K and AV « 6kms-i at Vlsr = 4.4 kms"' is de- 
tected. The CO 7 = 3-2 line is strong, 7mb ~ 10 K at the same 
velocity. A small CO map around IRS 46 reveals only a red 
wing, most likely associated with the outflow from IRS 44. 

High-resolution L- and M-band spectra {R k, 25,000) were 
obtained with the NIRSPEC spectrometer at the Keck II tele- 
scope^ to provide kinematic information and confirmation of 
the Spitzer HCN detection through the CO = 1 - 4.7 ^m 
and HCN (^3 = 1-0 3.0 /xm stretching mode ro-vibrational 
bands (see Figure 1^. The '^CO absorption lines are un- 
resolved at the quiescent cloud velocity of Vlsr ~ 4km s"'. 
However, the resolved '^CO (AV « 30kms-') and HCN 
(Ay w 20km s"') absorption lines are observed to be shifted, 
at Vlsr ~ -20kms"'. Of these no counterpart is observed in 
the HCN 7 = 4- 3 and CO 7=3-2 JCMT spectra. We asso- 
ciate the HCN, C2H2, and CO2 absorptions seen in the Spitzer 
spectra with the blue-shifted components. 

3. ANALYSIS 

The Spitzer spectra are analyzed using a pure absorption 
model assuming local thermodynamic equilibrium (LTE) ex- 
citation of the levels at a single temperature. The adopted 
metho d is described in detail in Lahuis & van Dishoec^ 
I2OOO) and Boonman et al. (2003), which includes references 
to the molecular parameters and data used in the model. The 
main fit parameters are the average temperature and integrated 

^ The JCMT is operated by the JAC in Hilo, Hawaii on behalf of PPARC 
(UK), NRC (Canada) and NWO (Netherlands) 

' The W.M. Keck Observatory is operated as a scientific partnership among 
Caltech, the University of California and NASA, made possible by the W.M. 
Keck Foundation 
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column density along the line of sight for a given intrinsic 
line width, defined by the Doppler fe-value. The resolved 
HCN 1/3 lines can be the result of multiple unresolved compo- 
nents along the line of sight, so that the intrinsic b value can 
be smaller than the observed line width. For small Zj-values 
< 2 km s"' , no good fit can be made to the C2H2 1^5 and CO2 1^2 
profiles. Therefore, b is taken to range from 2 to 12 km s"' . 

The best fits to the €2^21^5, HCN 1^2, and C02J^2 bands 
observed in the IRS spectrum give of ~ 700, 400, and 
300 K and column densities of 3, 5, and 10 xlO'^cm"^ re- 
spectively for b « 5kms"'. The uncertainty in b results in 
an uncertainty of 25 % in these values. The blue-shifted 
CO 1/ = 1 - absorption band gives Tex = 400 ± 100 K and 

= (2 ± 1) X 10'** cm"- corresponding to a minimum H2 col- 
umn density of 1 x lO^^cm"^ assuming a CO abundance of 
2 X 10"'* (all gas-phase carbon in CO), as appropriate for 
warm dense gas. The 9.7 /xm silicate depth corresponds 
to A^H = 3 X 10^^ cm"^ assuming a conversion factor: A^h = 
T9 7 X 3.5 X 10^^ cm"^ (see Draine 2003il). X-ray observations 
give A^H = 11(±7) X 10^^ cm"^ dlmanishi. Kovama & Tsuboil 
12001'). Assuming most hydrogen is in H2 and allowing for 
some foreground absorption in the latter two determinations, 
all estimates are consistent with A^(H2)=1 x 10^^ cm"^ within 
a factor of two. The resulting abundance estimates are 3, 5, 
and 10 xlO"*" for C2H2, HCN, and CO2, respectively. The 
density of the gas is at least lO^cm"^, required to thermalize 
HCN and C2H2. 

Additional constraints can be obtained from a combined 
analysis of the HCN 1^2 and ly^ bands (Figure |3}- Assuming 
the same excitation temperature for both bands, the required 
column density to fit the 1^3 band is higher by a factor of 4 than 
that found from the 1^2 band. This suggests that geometrical 
effects and emission filling in the absorption may play an im- 
portant role (see Boonman et al. 2003). One possibility is that 
the absorbing region is smaller than the continuum emitting 
region, the size of which may depend on wavelength. Simi- 
lar increases may be expected for C2H2 and CO2. Thus, the 
above cited abundances are lower limits. The inferred HCN 
abundance, with respect to both H2 and CO, is up to four or- 
ders of magnitude larger than that found in cold interstellar 
clouds. 

4. DISCUSSION 

Where does this hot gas rich in organic molecules reside? 
The first clue comes from the HCN submillimeter JCMT 
spectrum. For a HCN column density of > 10'^ cm"^ derived 
from the IR data, the 7=4-3 pure rotational line is highly opti- 
cally thick. Thus, Tmb is expected to be close to the excitation 
temperature of ~ 400 K if the emission would fill the beam. 
Although a weak emission line is observed at Vlsr ~ 4km s"', 
the broad -20km s"' component is not detected with a 3a 
limit of 0.02 K in a 1 km s"' bin. This gives a beam dilution 
> 2 x lO'* which, for the JCMT beam size of 15", implies a 
source diameter for the hot gas of < 0.1 1" or 13 AU diameter 
at the distance of Ophiuchus. A similar limit follows from the 
lack of a blue wing on the CO 7=3-2 line. 

The second clue comes from the high temperatures and den- 
sities of the molecular gas. In general, temperatures of a few 
hundred K are found in YSO environments only in the inner- 
most part of envelopes or in the inner disks. The velocity of 
the hot gas provides a final clue. The radial velocity of IRS 46 
is unknown, and it is possible that IRS 46 itself is at a velocity 
of-20kms"' cDoppmann et al. 20 05). More likely, however, 
IRS 46 is close to the nominal cloud velocity of 4.4 kms"' and 
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Fig. 4. — 2-D distribution of tlie temperature (left) and density (right) in the 
inner disk for the best-fitting model to the SED. Included are the r = 1 pho- 
tosphere at 14 fim (dot-dashed curves) and two lines of sight (dashed lines) 
at the disk inclination of 75°. Most of the molecular absorption originates 
along the lines of sight toward the photosphere in the hot (> 300 K) inner rim 
and toward the hot surface (1500 K) of the inner rim on the opposite side of 
the star (see Sectionl4l. 

the hot gas is blue-shifted by ^ 25 kms"'. 

Based on these arguments, three possibilities are considered 
for the location of this hot gas: (i) the inner layer of any rem- 
nant collapse envelope on scales < 100 AU around IRS 46; 

(ii) the inner < lOAU regions of a nearly edge-on disk; and 

(iii) dense hot gas at the footpoint of a wind launched from 
the inner disk. The first two options are examined through 
radiative transfer models to constrain the physical parameters 
of the source environment. 

To explore the first option, a spherically symmetric model 
with a power-law density distribution was constructed repro- 
ducing the SCUBA submillimeter continuum map and the 
SED of IRS 46 from mid-IR to submillimeter wavelengths 
following J0rgensen, Schoier & van Dishoeck (2002). This 
model has two severe problems: to reproduce the mid-IR con- 
tinuum emission, the temperature in the innermost envelope 
cannot exceed roughly 100 K; and a ^ 0.01 Mq envelope (es- 
timated from the submillimeter continuum data) produces a 
significant scattering nebulosity at near IR wavelengths which 
is no t observed in VLT-ISAAC images JPontoppidan et all 
l2005h . The envelope model also cannot explain the line 
widths and velocities unless IRS 46 itself is at -20 km s"' . 
Thus, most of the IR and submillimeter continuum emission 
must arise from the disk around IRS 46. 

To investigate the second option, the physical structure of 
the IRS 46 disk was constrained from th e observed SED in 
a manner similar to CRBR 2422.8-3423 JPontoppidan et alJ 
2005) using the 2-D axisymmet ric Mo n te C arlo radiative 
transfer code of Dullemond & DominikI J2004I) . Figure |4] 
shows the temperature and density structure of the hot in- 
ner part of the nearly edge-on self-shadowed flaring disk and 
Figure n the best fitting SED. There is some degeneracy in 
the parameters of the best-fitting model (see the discussion 
in Pontoppidan et al. 2005), but spatially resolved data are 
needed to further constrain the fits. The fit to the silicate fea- 
ture is sensitive to the assumed opacities, emission from the 
outer disk, and the presence of foreground absorption; these 
uncertainties have little impact on the part of the disk model 
relevant to this work, i.e., the dense inner disk region. 

The main continuum contribution from the disk at 3 - 14 /im 
comes from the puffed up inner rim and inner rim wall on 
the far side of the star. In the disk model, the integrated col- 
umn of dense gas (> 10*^ cm"-') toward these areas ranges from 
1 - 2 X 10^^ cm"^ with average temperatures of 300- 1500 K. 
This is consistent with the H2 column density and tempera- 
tures derived from the observations (see Section |3}- Indeed, 
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more generally, the observed temperatures are consistent with 
disk models which explicitly calculate the gas temperature 
in the inner disk. The gas temperatures may be even higher 
than the dust temperatures in t he upper layers (Glassaold et ^ 
12004; Gorti & Hollenbach 2004). The velocity broadened CO 
and HCN profiles (Figure |3} could result from absorption in 
the Keplerian inner rim at 0.1 AU (Vkep — 70kms"' for a 
0.5 M© star). If IRS 46 is at the cloud velocity of +4 km s~\ 
the blue-shifted absorption implies a deviation from Keplerian 
rotation in the disk plane, for example as the result of a close 
binary. A near edge-on disk explanation is consistent with a 
detection toward only one in a hundred objects, since a very 
small fraction of sources should have the right orientation. 

A blue-shifted velocity may also indicate that the observed 
absorption features originate in a wind emanating from the 
inner disk. Possibilities include a magnetocentrifugal wind 
(either an X-wind launched within 0.1 AU or a disk wind 
launched further out), a photoevaporative flow or stellar wind 
interacting with t he upper layer an d entraining molecular ma- 
terial (see e.g. Eisloffel et al. 2000). For the gas to be seen 
in absorption it must be comparable in size to the inner disk 
region responsible for the 14/im background, i.e. a few AU 
in radius. For a smaller region the background continuum 
will dominate whereas for a much larger region line emission 
from the warm molecular gas will fill in or dominate the ab- 
sorption. For a few AU region, the mass loss rate would be 
of order 10"^- 10"^ Mq yr"' assuming a density of > 10** 
cm~^, a total H2 column of 10^^ cm"-, and a flow velocity of 
~ 25 km s"' . IRS 46 does not show strong accretion signatures 
so a higher flow rate seems unlikely. Quantitative predictions 
in terms of velocities, column densities, densities and temper- 
atures of the molecular gas are needed to distinguish between 
the above models. 

Regardless of the precise origin in the disk or disk wind, the 
high inferred excitation temperatures of 400 - 900 K and high 
abundances of HCN and C2H2 of 10"''- 10"^ are characteris- 
tic of high temperature chemistry. Hot chemistry in general 
is dominated by evaporation of the molecules from the grains 
with subsequent gas-phase processing. At high temperatures. 



the hydrocarbon and nitrogen chemistry are enhanced as most 
of the oxygen is converted into water by neutral-neutral reac- 
tions. The abundances of molecules such as C2H2, CH4, and 
HCN can be increased by or ders of magnitude (e.g. Dotv et al] 
'2002; Rodg ers & Cha rnlev'2003) while at the same time the 
formation of CO2 is reduced as its primary formation route 
through OH is blocked. Interestingly, CO2 has a lower excita- 
tion temperature in our observations. The most recent models 
of inner disk chemistry pre dict enhanced abundan ces of HCN 
and C2H2. In particular, Mar kwick et alJ J2002h give HCN, 
C2H2, and CO2 abundances of 10"'' ~ 10"^ in the inner 1 AU 
of a protoplanetary disk, in good agreement with the abun- 
dances found in this work. 

In summary, we present the first detection of gaseous 
molecular absorption bands with Spitzer toward a solar-mass 
YSO, which offer direct, unique probes of hot organic chem- 
istry in its immediate environment. In addition, they provide 
independent constraints on the gas temperatures and velocity 
patterns. The most plausible scenario is that the absorption 
originates from the inner few AU of the circumstellar disk, 
perhaps at the footpoint of a disk wind. Further work, both 
observationally and theoretically, is required to prove this, in- 
cluding a determination of the velocity of IRS 46 itself and 
monitoring of the infrared lines at high spectral resolution 
to check for time variability of the radial velocity and/or ab- 
sorbing column along the line of sight. This detection of- 
fers prospects for future high spectral and spatial resolution 
mid-infrared and submillimeter searches for these and other 
organic molecules in emission in more face-on disks. 
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